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When quark-gluon plasma emerges in the wake of a heavy-ion collision, magnetic field cre¬ 
ated by the valence charges has already permitted the entire interaction region. Evolution 
of this “initial” field in plasma is governed by Maxwell equations in electrically conducting 
medium. As the plasma expands external valence charges induce magnetic field that also 
contributes to the total magnetic field in plasma. We solve the initial value problem describ¬ 
ing these processes and argue that the initial magnetic field often dominates over the one 
induced by the valence charges. In particular, it grows approximately proportional to the 
collision energy, unlike the induced component, which is energy-independent. As a result, 
magnetic held has a signihcant phenomenological inhuence on quark-gluon plasma at the 
LHC energies over its entire lifetime. 

I. INTRODUCTION 

In this paper we revisit the problem of magnetic field created by electrical currents of colliding 
relativistic heavy ions m- Since these currents experience very little deflection in the course of 
collision mm, the corresponding magnetic field depends on energy and geometry of the collision, 
and implicitly on the strong interaction dynamics through the electrical conductivity of the quark- 
gluon plasma (QGP) [21 [6]. Another important aspect, which is the main focus of this study, is 
transition dynamics from magnetic field in vacuum to the one in medium. To begin with, assume 
that QGP forms instantly at time t = to, where t is counted from the collision time in the laboratory 
frame. This time emerges in phenomenological models of QGP that favor rather small values as 
compared to the perturbation theory expectations, see e.g. HU- The earliest possible value of to 
is determined by the saturation momentum Qg as 1 /Qg and represents the time it takes to release 
most particles from the ion’s wave functions. At RHIC ~ 0.2 fm. At t < to we are dealing with 
electromagnetic field created by the valence charges in vacuum. Its magnetic component is given 
by the well-known formula ([^. At time t = to, when the QGP emerges, magnetic field permits 
the entire plasma. Starting at t = to and on behavior of magnetic field is governed by the Maxwell 
equations in plasma. These equations describe evolution of magnetic field in electrically conducting 
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QGP starting from its initial value at t = to- This component of the total magnetic field is referred 
to below as the “initial” magnetic field Another contribution to magnetic field is induced 

by valence charges moving outside of QGP and is referred to below as the “valence” contribution 
In previous publications the role of the initial field has not been properly recognized. In this 
paper we fill this void and moreover, argue that in most cases the main contribution stems from 
the initial field. 

The paper is organized as follows. In Sec. [IlJ-Sec. W we deal with magnetic field produced by 
a single point charge. In Sec. |l|] we consider magnetic field in vacuum, while in later sections - 


in electrically conducting QGP. The main result is given by equations (34), (35) which represent 
contributions of valence charges and the initial field respectively. A more realistic geometry is 
considered in Sec. [V| where we discuss the case of two electric charges colliding at a given impact 
parameter b. We also discuss there the effect of time dependent electrical conductivity on the 


magnetic field evolution. We discuss the results and summarize in Sec. VI 


II. MAGNETIC FIELD IN VACUUM 

In a relativistic heavy-ion collision, electromagnetic field is created by Z electric charges of one 
ion and Z electric charges of another ion moving in the opposite directions along, say, z axis such 
that ion centers are at a distance b away. Due to the superposition principle, the total classical 
field is a sum of fields of all charges. Thus, in order to find the total field it is sufficient to solve 
for a single electric charge e. In this section we briefly review a textbook case of electromagnetic 
field created in vacuum by a uniformly moving point charge e. Our intent here is to introduce 
notations, definitions etc. 

Before the QGP formation, viz. at f < to, the vector potential A^(r, t) of a point charge e moving 
along the trajectory z = vt satisfies the following equation 

V‘^Ai{ic,t) = d^A^{T,t), ( 1 ) 

where the electromagnetic current density due to a valence charge e is 

1 = evz6{z — vt)6{h). (2) 

* To avoid confusion we emphasize that both components are ultimately related to electrical charges of heavy-ions. 
The distinction only concerns our treatment of magnetic field at t > to as will be explained in detail in the 
forthcoming sections. 
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The momentum space representation is defined as follows 


. / .\ f d k j]c-r■ f d k^dk^ .K i 

= J (2^^ "ik. = J • 

With this normalization the Fourier component of the current reads 

It follows from 0 that the vector potential generated by the current 0 is 

lirevz Inevz 

- ik j^2 _ j^2y2 k^/'y‘^ + k'j_ 


In the configuration space we obtain 


Ai(r,f) = 


jevz 


where 7 = (1 — The corresponding magnetic field 

’-yevcf) b 


Bi = -dbAicf) = 


(3) 


(4) 


(5) 


( 6 ) 


(7) 


47r (6^ + 

This solution is valid until t = to at which time existence of electrically conducting medium must 
be taken into account. 


III. EXACT SOLUTION FOR CONSTANT ELECTRICAL CONDUCTIVITY 

Maxwell equations can be solved exactly for t > to in the case of constant electrical conductivity 
a. The vector potential A 2 satisfies the following equation 

V 2 A 2 (r, t) = A 2 (r, t) + adtA^iv, t) - j(r, t) ( 8 ) 

with the initial conditions 

A 2 (r, to) = Ai(r, to) = i4>(r, to), (9) 

'9tA2(r,to) = 9tAi(r,to) = ^4'(r,to) • (10) 

We assumed that permittivity and permeability of QGP is trivial. One can take a more accurate 
account of medium properties, which would yield more elaborate initial conditions. However, they 
are not expected to significantly change the final result. 
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In momentum space Eq. Q and the corresponding initial conditions Q,(10) read 

= a?Aj];(t) + a8iA2ij(() - etize"’**”', (11) 

A,k{<0) = z-fkdo) = . (12) 

evz 


= z4'i^(to) = -ikzV- 


^-ikzvtfj 


kill‘d + k‘j_ 

To solve 0 , we first consider the corresponding homogeneous equation 

- k'^a-^{t) = dia]^{t) + adta]^{t) , 


(13) 


(14) 


Seeking its solution in the form (x e we find, upon substitution into (14), that oj must obey 
one the following dispersion relations 


^cr 


id = UJ± = —— ± \ l k^ -r- 


cr^ 


(15) 


Thus, the general solution of the homogeneous equation (14), which describes propagation of the 
initial conditions, reads 


xjj(t) = 


(16) 


where a and fi are constants to be determined from the initial conditions (12) and (13). The 
particular solution due to the external current density is of the form (x where 6 is 

found upon substitution into 0 : 

ev 


S = 


— ikzva 


Thus, the general solution to (11) is 






ev 


—ikzvt 


k‘^ — k'iv'^ — ikzVa 


(17) 


(18) 


Applying the initial conditions (12) and (13) we can fix a and (3. The final result is 

~ ~ kzV -ikzvto 


"^2k ^ s ^ 


+ - 


a;+ — W- 
^ A>i 


+ w 


z(w+ — U!-) - 

1 


—loj-e 


0 )+ — UJ- 


+ e 


—ikzVt 


-T,. 




+ e 


(t—to) 


(19) 


z(w+ — UJ-) - L 

Fourier transformation to the configuration space yields exact analytical solution to the initial 


value problem ([^-( 10). Analytical and numerical evaluations of the integral over k are challenging. 
Fortunately, in the ultra-relativistic limit 7 ^ 1, which is relevant for the relativistic heavy-ion 


collisions, expression for the vector potential (19) significantly simplifies [B] . This is the subject of 
the next section. 
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IV. DIFFUSION APPROXIMATION 


For an ultra-relativistic charge moving along the trajectory z = vt, ^ k'j_,akz 

implying that (fsl) can be approximated by 


ViA 2 (r, t) = adtA 2 iz, t) - j(r, t) 


This equation, being of the first order in time derivative, requires only one initial condition 


( 20 ) 


A2(r,to) = Ai(r,to) = 


( 21 ) 


We can solve the initial value problem (20)-(21) for an arbitrary time-dependence of the conduc¬ 


tivity cr(t). Introducing a new “time”-variable A according to 

dt' 


m = 


'^0 


a{t') 


( 22 ) 


and transferring (20) to the momentum space we obtain 


^-L^2k ~ ik • 


(23) 


The corresponding homogeneous equation (i.e. (23) with = 0) is solved by 


a]^(A) = , 


(24) 


where (7 is a constant. To derive a particular solution, we treat C as a function of A and plug into 
We get 


C = ev r + D . 

Jo 


(25) 


Substituting into (24) we find the general solution to (23) 

nA 


^2 




(26) 


Since A(to) = 0, the initial condition (21) implies that D = 4>j^(r,to)- So finally. 


^ 2 . 


]j(t) = z Ieue-^" ^ $j^e-fciA(t')\ 

I Jto '^(^ ) “ J 


(27) 


In a particular case of constant electrical conductivity (|27|) simplifies to 

ev 1 




^ — ikzV 


(^-ikzVt _ ^-’l^[t-to)^-ikzVto'^ ^ 




(28) 
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This expression can be derived directly from (19), bnt the approach described in this section is 


more straightforward. Fourier transformation to the configuration space 


A 2 (r,t) = 


( 2 vr) 


2tt 


-e --L - 




(29) 


can be done using the following integrals 


(fk± 

r+00 

dk^ , 


to 

c 

' —00 

2 tt 


cfk± 

r-l-oo 

dk, ^ 


(2vr)^ . 

' —00 

2 tt 


yen 7 

4vr Jo 

00 

dk_\_Jo{k±_h) 

e 


exp ■ 


b 2 


4 [A(t)—A(t 


')]} 


ev 




4[A(t) - A(f')] 

„-ikzVto 


6{z — Vt') , 


(30) 


(31) 


Substituting (27) into (29), doing integrals (30),(31) and then integrating over t' yields 


A 2 (r,t) = 


ze 


exp 


-} 


4cj(z/f) 4[A(t) — A(z/n)] 


9{tv — z)9{z — vto) 


Jo 

Magnetic field can be calculated as in ([^ with the following result 

82 8val T Sinit J 

where the “valence” B^g^j and “initial” components are given by 

aTrb f 6 ^ 


(32) 


(33) 


eBvai(r,t) =0 


2a{z / v)[\{t) — A(z/n)]' 


exp 


4[A(t) - X{z/v)] 


9{tv — z)9{z — vto ), (34) 

(35) 


eBinit(E) ^) =070^^ / dA;_i_A;_L Ji(fe_L 6 ) exp /i;^A(t) — A:_l 7 |z — ntol} • 

Jo 

The fine structure constant a = e^/( 47 r). Note that at t = to, Bvai vanishes whereas Binit yields 
the initial condition 0- Binit is the field that permits the plasma as it emerges at t = to (at which 


time it coincides with Bi) and spreads in it according to (35). Unlike BvaO it strongly dependences 
on the collision energy 27 (in units of proton mass). 

Bval describes induced electromagnetic field generated as a response of QGP to electromagnetic 
field of the valence charge and builds up starting from t = to- Because of the two step-functions in 


(34) that reflect causality, Bval is finite only in the interval vto < z < vt. In particular, it vanishes 


in the midrapidity z = 0. At fixed z satisfying z > vto, Bval emerges when t = z/v. An important 
property of Bval is that its magnitude is independent of energy (since 1 ). 
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At early times since the QGP creation viz. t > to, expression in the exponent of (35) is such that 


k‘j_X <C A:_l 7 |z — vto\ implying that However, at later times when k‘j_X S> A:_l 7 |z — vto\, 

we get 


eBinit = 0 


'-yavbyjT: 
8A3/2 " 


e 8 a 




8A/ 


V8A 


(36) 


Since k±b ~ \/8 (which can be seen from Ji series expansion) and A ~ (t — to)/a we estimate that 


(36) is valid at times t satisfying 


^ “ ^0 1 ,, 

1-TT ^ - 1 = 1 ^ 0 . 

\z-vto\ 


(37) 


At z = 0, 6 = 7 fm and to = 0.2 fm this implies t ^ 1 fm, where I used a = 5.8 MeV known 


from the lattice calculations m, see also [I1HI6]. Furthermore, since 6^/8A <C 1 we expand ( |36[ ) 
to obtain the late-time behavior of the initial magnetic field 

jav,/7rb 2 


eB 


- init 


8A3/2 


(38) 


For constant a the late-time dependence (viz. t ^ to) is Hinit ~ l/t3/2. Notice that at late times 
the “valence” contribution decays as ~ 1/f^. It therefore emerges that the initial magnetic 
field dominates at early and late times. 

V. MAGNETIC FIELD OF TWO COUNTER-PROPAGATING CHARGES 


To calculate magnetic field in a heavy-ion collision one considers two sets of Z counter- 
propagating electric charges distributed according to one of the known nuclear density parameter- 
izations, see e.g. [3]. However, to study the time evolution of magnetic field it suffices to consider 
just two counter-propagating charges. The geometric symmetry of this configuration is similar to 
that of the event-average over many heavy-ion collisions at impact parameter b, but drastically 
reduces the computational time. The configuration that we consider is depicted in Fig. 


Let HB)(rj^,t) and B'^‘^\T 2 ,t) be magnitudes of the fields of the two charges, each given by (33) 


(35). We can express coordinates of the observation point relative to each charge r^ = b^ -|- zzi 


and 12 = b 2 -|- zz 2 in terms of their center-of-mass in cylindrical coordinates R, z, (j) as follows, see 

Fig.[T] 

72 sin (h 


ba = \/b^jA: -1-722-1- (—1)“ bRcoscp, tan (jia = 


72cos(/> — (—1)“ b/2 


Za = vt + (-1)'^ z . (39) 
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FIG. 1: Two counter-propagating charges e. One charges moves along the positive z-axis at z = vt, 
X = —6/2, y = 0 while another one moves in the opposite direction at z = —vt, x = 6/2, y = 0. 

where a = 1, 2 labels the charges. Noting that 6*-“^ oc (f)a and expressing (f)a in terms of b and (f) 
we obtain magnetic field in terms of the center-of-mass frame coordinates 

B sm{4> - (j)i) -|- sm{(j) - (1)2)] 

cos{(l) - (/i) -h B^‘^'>{r 2 ,t) cos((/> - (/ 2 )] , (40) 


where r^ and (j)a are replaced as indicated in (39). The result is shown in Fig. [^Fig. [^in terms of a 
dimensionless and unit-independent quantity eB/m^. In all figures impact parameter is 6 = 1 fm, 
observation point is at </> = 7 r/ 2 , i? = 7 fm (i.e. x = 0 and y = 7 fm), and 7 = 100 (except Fig. Q. 
Also indicated is the pseudorapidity 7 = — \tl[—{z/R) + \/ (z/i?)^ -|- 1]. Solid lines indicate the total 
magnetic field B, dashed lines represent the contribution of the initial condition Sinit and dotted 
lines stand for the contribution of the valence charges Sval- As discussed at the end of the previous 
section valence charge contribution decreases with time faster than that of the initial condition. 


eBIm^ eBIm^ 



FIG. 2: Magnetic field in units of vn?^le. a = 5.8 MeV, z = 0 fm [ij = 0). Left panel: to = 0.2 fm, right 
panel: to = 0.5 fm. Valence current does not contribute at all (Svai = 0). 
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eBIm^ eBIm^ 



FIG. 3: Magnetic field in units of m^je. a = 5.8 MeV, z = 0.6 fm (ry = 0.086). Left panel: to = 0.2 fm, 
right panel: to = 0.5 fm. Solid, dashed and dotted lines stand for B, Binit and Bvai- 

Fig. [^Fig.[^ depict magnetic field at constant electrical conductivity u = 5.8 MeV [l3]. In Fig.[^ 
we compare magnetic field that is generated when QGP emerges at to = 0.2 fm and at to = 0.5 fm. 
Since magnetic held in vacuum decreases as 1/t^, see Q, the late emergence of conducting medium 
means that the magnitude of the held in the former case is about 15 times larger than in the later. 
In both cases time-dependence of magnetic held in plasma is mild. Because of the step functions 
in pd] ) magnetic held at midrapidity z = 0 is entirely due to the initial held Binit- 

Fig. [3] is similar to Fig. [^except that z = 0.6 fm unlocking the “valence” contribution. Being 
independent of the initial value of magnetic held at to the “valence” contribution rapidly increases 
to its maximal value, that can be determined from ( |34[ ) |12j . It then decreases at larger t and 
becomes smaller than Binit- Sharp lines seen in Fig. [^indicate that the transition dynamics near 
t = to is not fully captured by the diffusion approximation. 

Energy dependence of magnetic held between the RHIC and LHC energies can be seen in 
Fig.i Binit grows approximately proportional to the collision energy 7 , whereas Bvai is energy 
independent. Thus, at the LHC magnetic held induced by valence charges is negligible. 

So far we considered only the case of constant electrical conductivity. In practice, however 
electrical conductivity is time-dependent. To see the impact of a time-dependence on the time 
evolution of magnetic held we consider two models. In model A we assume that QCP emerges 
instantly at t = to with cr = 5.8 MeV and then cools down as it expands according to the Bjorken 
scenario m- Namely, expansion is supposed to be isentropic nV = const, where n is the particle 
number density and V is plasma volume. Since n ^ and at early times expansion is one¬ 
dimensional V ~ t it follows that T oc Since a{t) oa T we conclude that a{t) ~ Thus 
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eBIm^ 


eBIm^ 




FIG. 4: Magnetic field in units of mj/e. a = 5.8 MeV, z = 0.2 fm to = 0.2 fm. Solid, dashed and dotted 
lines stand for B, Binit and i?vai- Left panel: 7 = 100 (RHIC), right panel: 7 = 2000 (LHC). 


a reasonable model for time dependence of electrical conductivity is 

= 2-V3(l I t/(„)./3 • ““'I’*' («) 

Another possibility is that the QGP does not appear as a thermal medium right away at t = to, 
rather it takes time r until the conductivity reaches its equilibrium value a. This can be described 
as 


a{t) =a(l-e , 


Model B. 


(42) 


We set conservatively r = 1 fm. Note that we cannot let cr{t) vanish at t = to because that would 


violate the diffusion approximation that lead to (20). However, (42) insures that cr(fo) ^ <7. 


In Fig. we contrast the two models. Similar calculation at constant conductivity is shown 
in the left panel of Fig. We observe that time-dependence (model A) significantly reduces 


magnetic field at later times. As far as model B is concerned, time dependence (42) affects mostly 


Hvai because it directly depends on cr{t), whereas Binit depends only on X{t), see (34),(35). Model 
B has minor effect on the total magnetic field, although one can certainly find regions in space-time 
where its effect is more pronounced. What actually matters is the initial time to at which one can 
treat the produced particle system as a medium. As long as conductivity is large enough at later 
times, magnetic field is fairly insensitive to the precise QGP dynamics. 


VI. SUMMARY 

Just before the QGP emerges, the interaction region is permitted by the primordial electromag¬ 
netic field created by valence charges of two heavy-ions. At the initial time to this magnetic field 
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eBIm^ eBIm^ 




FIG. 5: Magnetic field in units of rri^le. z = 0.2 fm to = 0.2 fm. Left panel: model A. Right panel: model 
B. Solid, dashed and dotted lines stand for B, Binit and B^ni- 


smoothly connects to the magnetic field in plasma and evolves according to the Maxwell equations 
in the electrically conducting medium. In addition to this “initial” magnetic field, there is another 
“valence” contribution that arises from the external valence electric charges inducing currents in 
the QGP. It has been tacitly assumed that the former contribution is not important [^. In this 
paper we argued to the contrary, that the initial magnetic field dominates at very early and later 
times and increases much faster with the collision energy than the “valence” contribution. 

We also studied the effect of time dependence of electrical conductivity and concluded that at 
early times it has a rather minor effect on the field strength, as long as the produced particle system 
can be treated as a medium at early enough time. However, towards the later times of plasma 
evolution, time-dependence of electrical conductivity plays an important role. In the Bjorken 
scenario it leads to much weaker fields as compared to the constant conductivity case. 

We considered the case of two counter-propagating charges that gives an accurate picture for 
the time dependence of the event-averaged fields in heavy-ion collisions. Scaling the result with 
Z we can obtain an estimate of the magnetic field strength in heavy-ion collisions. Calculating 
the spatial distribution requires an accurate account of the exact nuclear geometry, which is not 
difficult using the results reported in this paper. 
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